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Abstract 

A self-action problem for a pointlike charged particle arbitrarily moving in flat spacetime 
of three dimensions is considered. Outgoing waves carry energy-momentum and angular 
momentum; the radiation removes energy and angular momentum from the source which 
then undergoes a radiation reaction. We decompose Noether quantities carried by electro- 
magnetic field into bound and radiative components. The bound terms are absorbed by 
individual particle's characteristics within the renormalization procedure. Radiative terms 
together with already renormalized 3-momentum and angular momentum of pointlike charge 
constitute the total conserved quantities of our particle plus field system. Their differential 
consequences yield the effective equation of motion of radiating charge in an external elec- 
tromagnetic field. In this integrodifferential equation the radiation reaction is determined 
by Lorentz force of pointlike charge acting upon itself plus nonlocal term which provides 
finitcncss of the self-action. 

PACS numbers: 03.50.De, ll.10.Gh, 11.30.Cp, ll.10.Kk 

1 Introduction 

Recently [1, 2], there has been considerable interest in renormalization procedure in classical 
electrodynamics of a point particle moving in flat space-time of arbitrary dimensions. The main 
task is to derive the analog of the well-known Lorentz-Dirac equation [3]. The Lorentz-Dirac 
equation is an equation of motion for a charged particle under the influence of an external force 
as well as its own electromagnetic field. (For a modern review see Refs. [5] and [6].) 

A special attention in Refs. [1] and [2] is devoted to the mass renormalization in 2 + 1 theory. 
(Note that electrodynamics in Minkowski space M3 is quite different from the conventional 
3 + 1 electrodynamics where one space dimension is reduced because of symmetry of a specific 
problem. For example, small charged balls on a plane are interacted inversely with the square 
of the distance between them, while in M3 the Coulomb field of a small static charged disk 
scales as |r| _1 .) An essential feature of 2 + 1 electrodynamics is that Huygens principle does 
not hold and radiation develops a tail, as it is in curved space-time of four dimensions [7] where 



'Electronic mail: yar@ph.icmp.lviv.ua 



1 



electromagnetic waves propagate not just at the speed of light but all speeds smaller than or 
equal to it. 

In Refs. [1] and [2] the self-force on a pointlike particle is calculated from the local fields in 
the immediate vicinity of its trajectory. The schemes involve some prescriptions for subtracting 
away the infinite contributions to the force due to the singular nature of the field on the particle's 
world line. In Ref. [2] the procedure of regularization is based on the methods of functional 
analysis which are applied to the Tailor expansion of the retarded Green's function. The authors 
derive the covariant analog of the Lorentz-Dirac equation which is something other than that 
obtained in Ref. [1]. Both the divergent self-energy absorbed by the "bare" mass of pointlikc 
charge and the radiative term which leads an independent existence are nonlocal. (They depend 
not only on the current state of motion of the particle but also on its past history.) 

In this paper we develop a consistent regularization procedure which exploits the symmetry 
properties of 2 + 1 electrodynamics. It can be summarized as a simple rule which obeys the 
spirit of the Dirac scheme of decomposition of the vector potential of a pointlike charge. 

According to the scheme proposed by Dirac in his classical paper [3], one can decom- 
pose retarded Green's function associated with the four-dimensional Maxwell field equation 
G ret (x, z) = G sym (x, z) + G rad (x, z). The first term, G sym (x, z), is one-half sum of the retarded 
and the advanced Green's functions; it is just singular as G ret (x, z). The second one, G rad (x, z), 
is one-half of the retarded minus one-half of the advanced Green's functions; it satisfies the 
homogeneous wave equation. Convolving the source with the Green's functions G sym (x,z) and 
G rad (x, z) yields the singular and the radiative parts of the vector potential of a pointlike charge, 
respectively. 

The analogous decomposition of Green's function in curved space-time is much more delicate 
because of richer causal structure. Detweiler and Whiting [9] modified the singular Green's 
function by means of two-point function v(x,z) which is symmetric in its arguments. It is 
constructed from the solutions of the homogeneous wave equation in such a way that a new 
symmetric Green's function G s (x,z) = G sym (x,z) + l/(8ir)v(x, z) has no support within the 
null cone. 

It is obvious that the physically relevant solution to the wave equation is the retarded so- 
lution. In Ref. [10] the Lorentz-Dirac equation is derived within the framework of retarded 
causality. Teitelboim substituted the retarded Lienard-Wiechert fields in the electromagnetic 
field's stress-energy tensor and computed the flow of energy-momentum which flows across a 
tilted hyperplane which is orthogonal to particle's four-velocity at the instant of observation. 
The effective equation of motion is obtained in Ref. [10] via consideration of energy-momentum 
conservation. Similarly, Lopez and Villarroel [11] found the total angular momentum carried by 
the electromagnetic field of a pointlike charge. It was shown [12] that the Lorentz-Dirac equation 
can be derived from the energy-momentum and angular momentum balance equations. In Ref. 
[13] the analog of the Lorentz-Dirac equation in six dimensions is obtained via analysis of 21 
conserved quantities which correspond to the symmetry of an isolated point particle coupled 
with an electromagnetic field. (First this equation was obtained by Kosyakov in Ref. [14] via 
the consideration of energy-momentum conservation. An alternative derivation was produced 
by Kazinski et al in Ref. [2].) 

In nonlocal theories, the computation of Noether quantities is highly nontrivial. Quinn and 
Wald [15] studied the energy-momentum conservation for a point charge moving in curved space- 
time. The Stokes theorem is applied to the integral of flux of electromagnetic energy over the 
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compact region V(t + ,t~). It is expanded to the limits t~ — > — oo and i + — > +00, so that finally 
the boundary of the integration domain involves smooth spacelike hypersurfaces at the remote 
past and in the distant future. The space-time is asymptotically flat here. The authors proved 
that the net energy radiated to infinity is equal to the total work done on the particle by the 
electromagnetic self-force. (DeWitt-Brehme [7, 8] radiation-reaction force is meant.) It was 
shown also [15] that the total work done by the gravitational self-force is equal to the energy 
radiated (in gravitational waves) by the particle. (The effective equation of motion of a point 
mass undergoing radiation reaction is obtained in Ref. [16]; see also review [18] where the 
motion of a point electric charge, a point scalar charge, and a point mass in curved space-time 
is considered in details.) 

In Ref. [17] Quinn derived the effective equation of motion of a point particle coupled with 
a scalar field moving in curved space-time. The author establishes that the total work done by 
the scalar self-force matches the amount of energy radiated away by the particle. 

In the present paper we calculate the total flows of energy-momentum and angular momen- 
tum of the retarded field which flow across a hyperplane S t = {y 6 M3 : y° = t} associated with 
an unmoving observer. The field is generated by a pointlike charge arbitrarily moving in flat 
Minkowski space M3 of three dimensions. This paper is organized as follows. In Section 2 we 
recall the retarded and the advanced Green's functions associated with the three-dimensional 
D'Alembert operator. Convolving them with the point source, we derive the retarded and the 
advanced vector potential and field strengths. In Section 4 we trace a series of stages in the 
calculation of a surface integral which gives the energy-momentum carried by the retarded elec- 
tromagnetic field. (Details are given in the appendixes.) We introduce appropriate coordinate 
system centered on an accelerated world line and we express the components of the Maxwell 
energy-momentum tensor density in terms of these curvilinear coordinates. We integrate it over 
the variables which parametrize the surface of integration £j . The resulting expression becomes 
a combination of two-point functions depending on the state of the particle's motion at instants 
t\ and ti before the observation instant t. They are integrated over the particle's world line 
twice. We arrange them in Section 3. We split the momentum three-vector carried by electro- 
magnetic field into singular and radiative parts by means of the Dirac scheme which deals with 
fields defined on the world line only. All diverging quantities have disappeared into the proce- 
dure of mass renormalization while radiative terms lead independent existence. In an analogous 
way we analyze the angular momentum of the electromagnetic field. Total energy-momentum 
and total angular momentum of our particle plus field system depend on already renormalized 
particle's individual characteristics and radiative parts of "electromagnetic" Noether quantities. 
Having differentiated the conserved quantities we derive the effective equation of motion of a 
radiating charge. In Section 5, we discuss the result and its implications. 



2 Electromagnetic potential and electromagnetic field 
in 2+1 theory 

We consider an electromagnetic field produced by a particle with 5-shaped distribution of the 
electric charge e moving on a world line (CM3 described by functions z a (r) of proper time r. 
The Maxwell equation 

F°<% = 2nj a (2.1) 
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where current density j a is given by 

/+oo 
dTu a (T)5^(y-z(T)) (2.2) 
-oo 

governs the propagation of the electromagnetic field. u a {r) denotes the (normalized) three- 
velocity vector dz a {r)/dT and 5^\y — z) = 5(y° — z°)5(y 1 — z 1 )5(y 2 — z 2 ) is a three-dimensional 
Dirac distribution supported on the particle's world line £. Both the strength tensor F a/3 and 
the current density j a are evaluated at a field point y £ M3. (We choose Minkowski metric 
tensor r] a n = diag(— 1, 1, 1) which we shall use to raise and lower indices. Greek indices run from 
to 3, and Latin indices from 1 to 2; summation over repeated indices understood throughout 
the paper.) 

We express the electromagnetic field in terms of a vector potential, F = dA. We impose the 
Lorentz gauge A a ^ a = 0; then the Maxwell field equation (2.1) becomes 

UA a = -2irj a . (2.3) 

In 2+1 theory the retarded Green's function associated with the D'Alembert operator □ := 
rj a ^d a dfs has the form [1, 2] 

Gg 1 fo,*) = 9(/ 7°- |y ; x|) . (2.4) 

9(y° — x° — |y — x|) is step function defined to be one if y° — x° > |y — x|, and defined to be 
zero otherwise. Synge's world function a(y, x) is numerically equal to half the squared distance 
between y and x: 

a(y,x) = ^r ]a p(y a -x a )(yP-xP). (2.5) 

The first is y, to which we refer the "field point", while the second is x, to which we refer the 
"emission point" . 

While in four-dimensional Minkowski space-time the retarded Green's function has support 
on the future light cone of the emission point x, in 2 + 1 electrodynamics its support extends 
inside the light cone as well. 

Using the retarded Green function (2.4) and the charge-current density (2.2) we construct 
the retarded Lienard-Wiechert potential in three dimensions. Denoting = y M — z M (r) as the 
unique timelike (or null) geodesic connecting a field point y to the emission point z(t) G C> we 
arrive at 

^« feHe £^o_ |K|) _^ m 

where the dot denotes the scalar product of three-vector K on itself (it is equal to double Synge's 
function of field point y and emission point z(t)). 

We now turn to the calculation of electromagnetic field FJ^ = d fJi A r u et — d v A r ^ 1 generated by 
an arbitrarily moving pointlike charge. It consists of two quite different terms. The first term 
is due to differentiation of ^-function involved in the vector potential (2.6): 

F$)= lim g u,K -u v K r 
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Figure 1: In four dimensions the retarded (advanced) field at observation point P(y) is generated by a 
single event in space-time: the intersection of the world line and P's past (future) light cone. In three 
dimensions the retarded field depends also on the particle's history before T ret (y). The advanced field 
depends on the particle's history after r adv (y) . The vector K is a vector pointing from the emission point 
z(t) € C to field point P. 



T ret (y) denotes the proper-time parameter at the point on the world line which links with y by the 
unique future-directed null geodesic. Since T ret (y) is the root of algebraic equation K° — |K| = 
the 5-term (2.7) diverges. 
The second term is 

= - e f +0 ° dT6(K° - \K\) UflKu ~ UuK>1 

r ret (y) v K - v K 

[-(K-K)]W K ! 

We see that the strength tensor FT^ of the adjunct electromagnetic field consists of terms 
proportional to 5— and functions: F ret = + p( e \ The terms separately are singular. 
The singularity, however, can be removed from the sum of F^ and p( e \ Using the identity 



[-{K-K)\W -(K-u) d Ty /-(K-K) 



(2.9) 



in eq.(2.8) yields 



F {6) e UjxKv_ ~ UyKji 

V-(K-K) -(K-u) 



(2.10) 
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rr ret (y) dr (a^Ky-auK^ u^K u - u u K^ \ 

+ e L ^qnyl -(*-,.) + \-(k-uw ' 1 + « A -"| 

after integration by parts. Summing up (2.7) and (2.10) and taking into account that 1 / — (K ■ K) 
vanishes whenever r —> — oo 1 , we finally obtain the expression 

which is regular on the light cone. (It diverges on the particle's trajectory only.) The symbol A 
denotes the wedge product. The invariant quantity 

r = -(K-u) (2.12) 

is an affine parameter on the time-like (null) geodesic that links y to z(t); it can be loosely 
interpreted as the time delay between y and z(t) as measured by an observer moving with the 
particle. When r = r ret (y), the parameter r is also the spatial distance between z(T ret ) and y 
as measured in this momentarily comoving Lorentz frame. 

In three dimensions the advanced Green's function 

o».^ 9( - / ;:°" |y r'" ) <»■«) 

y / -2o-(y,x) 

is nonzero in the past of x. The advanced strength tensor 

vw,,, \ f + °° dr (aAK u A K r , rr N1 ] 

F adv (y) = el - = = \— — + — —[l + (K-a)}\ (2.14) 



^{y) \J ~{K ■ K) 



is generated by the point charge during its entire future history following the advanced time 
associated with y (see figure 1). 



3 Equation of motion of radiating charge 

In this section we derive the "three-dimensional" analog of the Lorentz-Dirac equation via analy- 
sis of the energy-momentum and angular momentum balance equations. The momentum three- 
vector carried by the electromagnetic field is calculated in the next Section and in Appendixes 
D, E, and F; the total angular momentum is obtained in Appendix G. We split the Noether 
quantities into bound (singular) and radiative (regular) parts. The energy-momentum and 
angular momentum of bare sources absorb the bound terms within regularization procedure. 
Already renormalized characteristics of charged particles are proclaimed to be finite. Together 
with radiative terms, they constitute the total energy-momentum and angular momentum of 
our particle plus field system which are properly conserved. 

To find out electromagnetic field's energy-momentum, we integrate the Maxwell energy- 
momentum tensor density over the plane Y, t = {y £ M3 : y° = t}. The resulting expressions 

x We assume that average velocities are not large enough to initiate particle creation and annihilation, so that 
"space contribution" |K| can not match with an extremely large zeroth component K . 
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(4.3.6) and (4.3.7), presented in the next Section, can be rewritten in manifestly covariant 
fashion: 

AW = f f *»-#^=f (3.1) 



-oo 

fTl 



+ 



e 2 / d,T\ I d.T2 

J —oc J — 

f *> f 

J —OO J TO 



dTl 

OO J T 2 



(fa q 1 ' 1 d f \ 1 d i (< 



ari^ra (2a) 3 / 2 2 dn Vv 7 ^/ 2 9r 2 Vv 7 ^ 



(We omit structureless terms which arise due to choice of non-covariant surface of integration.) 
Index 1 indicates that particle's velocity or position is referred to the instant t\ g] — 00, r] 
while index 2 says that the particle's characteristics are evaluated at instant t 2 < t\. Here 
= — Z2 defines the unique timelike geodesic connecting a field point z(ti) £ ( to an 
emission point z{t 2 ) G £; by a we mean the Synge's world function (2.5) of z\ and z 2 , taken 
with opposite sign: 

v(Ti,T 2 ) = -±(q-q). (3.2) 

Two double integrals over (proper) time variables (one about the other) describe integration 
over the domain D T = {(ti,t 2 ) G M 2 : t\ g] — oo,r],T2 < ti}. 

We have to decompose the expression eq.(3.1) into singular and regular parts. Following 
Ref.[2], we postulate that splitting should satisfy the conditions: 

• proper non-accelerating limit of singular and regular parts; 

• proper short-distance behavior of regular part; 

• Poincare invariance and reparametrization invariance. 

By "proper short-distance behavior" we mean the finiteness of integrand at the edge r 2 = t\ of 
the domain D T . 

So, we take one-half of the first term in between the square brackets under the double integral 
signs in eq.(3.1): 

1 d 2 a q^ 1 (u\ ■ u 2 )q^ 



2dT 1 dT 2 {2afl 2 2[-( g - 9 )]3/ 2 

and add the second term 

Id ( U2 \ _ 1 (u\ ■ q)u^ 

2dT- 1 W2a z ) = 2[-(q-q)} 3/2 ' 

It is convenient to rewrite the resulting expression as follows: 



(3.3) 
(3.4) 



„ 1 -v%q» + u%q a 
V ^ = 2 U1 >" [-(q-q)]W • (3 ' 5) 

We introduce the function 

GiU(Ti) = e Ui, Q J dT 2 r_ ( )|3/2 ^ 
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which is the convolution u^(ti)F^ of three- velocity and non-local part (2.8) of the retarded 
strength tensor evaluated at point z{t\) G C- It is intimately connected with the retarded field 
(2.11) generated at point z{t\) G ( by the portion of the world line that corresponds to the 
interval — oo < r 2 < T\\ 

G? et (T 1 ) = ^a»-eu 1 , a F^(T 1 ). (3.7) 

(It may be checked via integration by parts.) 

Next we take the remaining one-half of the first term and add the third term: 

1 (ui ■ u 2 )q t " 1 (u2 ■ g)^i 1 -ufq^ + u^q a 

2[-(q. g )]3/2 + 2 [_(g . g )]3/2 " [-(q . g)]3/2 • I • > 
We change the order of integration of this term over D r : 

e 2 T . P , -<g^ + <£° e 2 /" r -iff+itf 

2L dTi L d ™« i-o-w = ^l-J™»L " Ti i^iw (3 - 9) 

and interchange indices "first" and "second". Taking into account that q(T 2 ,Ti) = — q(ri,T 2 ), 
we finally obtain: 

-iL dTiui «L dT2 [-(»-, ■ (3 - 10) 

The integrand coincides with that under integral sign in the right-hand side of eq.(3.6) while the 
domain of inner integration corresponds to the interval T\ < t 2 < r. We introduce the function 

G adv( r l' T ) = e «l,a / <fra r /. ..Ma/9 • t 3 " 11 ) 



which is the convolution u^(ti)F^ of three- velocity and non-local part of the advanced strength 
tensor evaluated at point z{j\) G £. If the instant of observation r tends to +00, this function 
can be rewritten as 

G*LM) = y< " ™ ha FZ(n). (3.12) 

The second term is convolution of velocity with the advanced field (2.14) generated at point 
z(ti) S ( by the portion of the world line that corresponds to the interval [ti,+oo[. The 
relations (3.7) and (3.12) are symmetric upon future and past. 

We see that the double integral in eq.(3.1) can be expressed as one-half of G ret minus one-half 
of G a dv integrated over the world line (. For this reason we proclaim the expression 

K(r) = \f dn[GJU(n)-G^ v (n,r)] (3.13) 



= e 2 I dri I dr 2 



Tl 



— OO 



(u 1 -u 2 )q fl 1 (v,!- q)u^ | 1 (u 2 ■ q)u% 



[-(q-q)?' 2 2[-(g.g)]3/2 2 [-{q • q)}^ 



the radiative part of energy-momentum carried by the electromagnetic field. The situation is 
pictured in figure 2. 
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We evaluate the short-distance behavior of the expression under the double integral in 
eq.(3.13). Let t\ be fixed and t\ — t 2 := A be a small parameter. With a degree of accu- 
racy sufficient for our purposes 



V-(q • q) 



A 
A 



(3.14) 



u i ~ a \ ^ + ^1 



A ,,A 2 ^ 

v 

A 2 



Substituting these into integrand of the double integral of eq.(3.13) and passing to the limit 
A — > yields regular expression 



lim 

T2~ >Tl 



1 



-u%qV + vT 2 q a 1 



w 



-(q-q)] 3/2 



-(ax) <-^<- 



(3.15) 



Therefore the subscript "R" stands for "regular" as well as for "radiative" . 
Alternatively, choosing the linear superposition 



Ps 



dn[GJU(n) + 0^(71^)] 



(3.16) 



we restore the first term in the right-hand side of eq.(3.1). Indeed, having integrated the half-sum 
of the functions (3.6) and (3.11) over £, we obtain 



e 2 r 

Mr) = Y j 



T /J- 

dr 2 - ; 

00 y / 2a(Tl,T 2 ) 



+ 



T\=T2 



f 

J — C 



v / 2cr(ri,r 2 ) 



T2- 



dr2- r 

00 y / 2cr(r,T2) 



(3.17) 



Since this non-local term diverges, in eq.(3.16) the subscript "S" stands for "singular" as well 
as "symmetric". 

In the specific case of uniformly moving charge ^2a(r, t 2 ) = t — t 2 . Hence Pg(r) coincides 
with that obtained in Appendix A where rectilinear uniform motion is considered (see eq.(A.14)). 
Since the bracketed integrand in (3.13) vanishes if = const, nonaccelerating charge does not 
radiate. 

We therefore introduce the radiative part of energy-momentum and postulate that it, and 
it alone exerts a force on the particle. The singular part ps should be coupled with the particle's 
three-momentum, so that "dressed" charged particle would not undergo any additional radiation 
reaction. The already renormalized particle's individual three- momentum, say p pa rt) together 
with pn constitute the total energy-momentum of our particle plus field system: P = p par t + Pr- 
Since P does not change with time, its time derivative yields 



Ppart( T ) 



~Pr 



(3.18) 



j — c 



ds 



-u"q^ + Us q' 
[2a(r,s)] 3 / 2 



-u^q^ + Urq' 
[2<t(t, s)] 3 / 2 . 
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Figure 2: We call "retarded" the term (3.6) with integration over the portion of the world line before t\. 
We call "advanced" the term (3.11) with integration over the portion of the world line after t\. For an 
observer placed at point z(r\) <E C the regular part (3.13) of electromagnetic field momentum looks as the 
combination of incoming and outgoing radiation, and yet the retarded causality is not violated. We still 
consider the interference of outgoing waves presented at the observation instant r. The electromagnetic 
field carries information about the charge's past. 

(The overdot means the derivation with respect to proper time r.) Here index r indicates that 
the particle's velocity or position is referred to the observation instant r while index s indicates 
that the particle's characteristics are evaluated at instant s < r. 

Our next task is to derive an expression which explain how three-momentum of a "dressed" 
charged particle depends on its individual characteristics (velocity, position, mass etc.). We do 
not make any assumptions about the particle structure, its charge distribution and its size. We 
only assume that the particle 3-momentum p pa rt is finite. To find out the desired expression we 
analyze conserved quantities corresponding to the invariance of the theory under proper homo- 
geneous Lorentz transformations. The total angular momentum, say M, consists of particle's 
angular momentum z Ap part and radiative part of angular momentum carried by electromagnetic 



field: 




(3.19) 



(Singular part is absorbed by p par t-) The last term is calculated in Appendix G: 



p 2 rr rTi 



11 V V IX , fX V V LI \ 

Z \V\2 — Z \V\2 Z 2V2\ ~ Z 2P21) 



(3.20) 



J—oo J —CO 



where two-point function is given by eq.(3.5). 

Having differentiated eq.(3.19) and inserting eq.(3.18) we arrive at the equation 
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where the symbol A denotes the wedge product. We obtain the system of tree linear equations 
in three components of the particle's momentum. Its rank is equal to 2. Therefore, an arbitrary 
scalar function to(t) arises: 



^{r)=mu^{T) + 



uPis) - uP(t) 

ds , — . 

oo \/2cr(r, s) 



(3.22) 



(We choose the simplest expression that is finite near the point of observation.) We see that, 
apart from the usual velocity term, the particle's 3-momentum contains also nonlocal contribu- 
tion from the particle's electromagnetic field. 

The scalar product of the particle three-velocity on the first-order time-derivative of the 
particle three- momentum (3.18) is as follows: 

o2 



(Ppart ■ U T ) = — 



ds 



(u T ■ u. 



(ut ■ q) (us ■ q) 
[2cr] 3 /2 + [ 2(J ] 3/2 



(3.23) 



Since (u ■ a) = 0, the scalar product of particle acceleration on the particle three-momentum 
(3.22) is given by 



(PI 



part 



ds 



[a T ■ u s ) 



2a 



Summing up eqs.(3.23) and (3.24) we obtain 



d_ 



part 



U-r 



ds 



d_ 



(ll T ■ U s ) 



2a 



(u s ■ q) \ 
[2a] 3 / 2 J " 



Alternatively, the scalar product of 3-momentum (3.22) and 3-velocity is as follows: 



(Ppart • Ut) = ~m + 



ds 



[u T ■u s ) + l 



(3.24) 



(3.25) 



(3.26) 



Further we compare its differential consequence with eq.(3.25). A surprising feature of the 
already renormalized dynamical mass m is that it depends on r: 



m 



-I is 



(q ■ u T ) - (q ■ u s ) 
[2a]3/2 



(3.27) 



It is interesting that a similar phenomenon occurs in the theory which describes a pointlikc 
charge coupled with massless scalar field in flat space-time of three dimensions [19]. The charge 
loses its mass through the emission of monopole radiation. 

Having integrated derivative (3.27) over the world line (, we obtain 



m 



to 



= mo + 



dr\ 



Tl 



dro 



d_ 

3t\ 



1 



d_ 

dr 2 



1 



(3.28) 



ds 



AM 



where too is an infinite bare mass of the particle. Inserting this into eq.(3.22), we arrive at the 
equality Pp art ( T ) = m oUr+Ps which shows that the particle's momentum renormalization agrees 
with the renormalization of mass. 
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The main goal of the present paper is to compute the effective equation of motion of radiating 
charge in 2 + 1 dimensions. To do it we replace Pp art in the left-hand side of eq.(3.18) by 
differential consequence of eq.(3.22) where the right-hand side of eq.(3.27) substitutes for m. At 
the end of a straightforward calculations, we obtain 



mai. 



F 

J — oc 



ds 



e 



2 



-u*q» + u>£q a 1 a% 
_ _ _ 



-a 



[2d] 3 / 2 
ds 



(3.29) 



The first term in the right-hand side of this equation looks horribly irrelevant. Relation (3.7) 
prompts that the retarded Lorentz self-force should be substituted for the combination of the 
first (local) and of the second (non-local) terms in the right-hand side of eq.(3.29). If an external 
electromagnetic field F ext is applied, the equation of motion of radiating charge in 2 + 1 theory 
becomes 

e 2 r dn 

ma» = eu T , a F^(r) + -< / / __ + eu T , a F^ t (3.30) 



where 



r ds f 

i-oc a/2ct(t, s) I 



-f 

J — I 



dsr a ( T ,s) 



2 I 1 + (°« " 9)1 + 



(3.31) 



is the field strengths at point z(t) G C generated by a portion of the world line before the 
observation instant r. The non-local term in eq.(3.30) which is proportional to the particle's 
acceleration a(r) arises also in Ref.[2]. It provides proper short-distance behavior of the radia- 
tion backreaction. If s — > r, the integrand tends to three-dimensional analog of the Abraham 
radiation reaction vector: 



lim 



eu T , a r a (r,s) + 



'-e 2 (d" - aV) . 



(3.32) 



(All quantities on the right-hand side refer to the instant of observation r.) 

If one moves the second term to the left-hand side of eq.(3.30), they restore unphysical 
motion equation which follows from variational principle: it involves an infinite bare mass and 
divergent Lorentz self- force. 

4 Energy- momentum of electromagnetic field in 2+1 dimen- 
sions 



In this section we trace a series of stages in calculation of the surface integral 



(4.1) 
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which gives the energy-momentum carried by electromagnetic field of a point-like source arbi- 
trarily moving in M3. In Appendix D, Appendix E and Appendix F we perform the computation 
in detail. 

In eq.(4.1) da^ is the vectorial surface element on an arbitrary space-like surface S. The 
electromagnetic field's stress-energy tensor T has the components 

2vrT^ = F ttX F v x - 1/A V ^F kX F kX (4.2) 
where F is the nonlocal strength tensor (2.11). 



4.1. Coordinate system 

In general, the rate of radiation does not depend on the shape of S. We choose the simplest 
plane S t = {y G M3 : y° = t} associated with unmoving inertial observer. If parametrization 
of the world line is provided by a disjoint union of planes particle's velocity takes the form 
= jv^jV^ = (l,z l ), and acceleration = 7 4 (vv)u' 1 + j 2 ^; factor 7 = — v 2 . (The 

overdot indicates differentiation with respect to t.) Electromagnetic field (2.11) takes the form 

r-ret, n f tret(y) dt fvf\K VAK, 2 fryr 

F iy)=e L 7^m\— + — [l +( K v)] \ (4L1) 

where = (0,v l ) and r = —{K ■ v). (Although we use the same notation, r should not be 
confused with manifestly covariant parameter (2.12).) 

Huygens principle does not hold in three dimensions and radiation develops a tail (see figure 
3). In 3D the circle C(z(0),t) = {y G M 3 : (y ) 2 = £ 2 =1 (y* - ^(0)) 2 ,y° = t} is filled up 
by electromagnetic radiation even if interval At — > 0. (The period of time during which the 
point source emanates is meant.) So, a point z(t±) G C produces the disk of radius t — t\ in 
the observation plane S< = {y G M3 : y° = t}. This property reflects the fact that in M3 the 
electromagnetic field at y is generated by the portion of the world line that corresponds to the 
interval — 00 < r < r re *(y); this represents the past history of the particle. 

We introduce coordinate system associated with two points on a particle's world line labelled 
by instants t\ and ti (see figure 4). Flat space-time M3 becomes a disjoint union of planes 
T, t = {y G M3 : y° = t}. A plane H t is a union of (retarded) disks centered on a world line of 
the particle. The disk 



C(z(t a ),t - t a ) = {y G M 3 : y° - t a > ~ ^(*a)) 2 , y° = t} (4.1.2) 

is bounded by the intersection of the future light cone generated by null rays emanating from 
z(t a ) G C m au possible directions, and plane T,f. The circular spot (4.1.2) is filled up by circles 
of radii R G [0, t — t a ] centered on points on a line connecting points z l {t\) and z l {t2)- Points in 
an R— circle are distinguished by polar angle (p. We define the coordinate transformation locally 
written as 

y° = t (4.1.3) 

y l = az i {t 1 )+pz i (t 2 ) + RuJ i j n j 
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Figure 3: Let the point source radiates within the interval [0, At]. In four dimensions the support of 
the Maxwell energy-momentum tensor density in hypcrplane y — t is in between two spheres centered 
at points z l (0) (cross symbol) and z l (At) (box symbol) with radii t and t — At, respectively. In three 
dimensions the radiation fills the disk with radius t centered at point z l (0) (cross symbol) even if the 
interval shrinks to zero. 



where a + (5 = 1 and n 3 = (cos if, sin ip). Orthogonal matrix uj is given by eq.(B.2) (see Appendix 
B). It rotates space axes till new y 1 — axis be directed along two-vector q := z(ti) — zfo)- 

The integration of energy and momentum densities over two-dimensional plane y° = const 
means the study of interference of outgoing electromagnetic waves emitted by different points 
on particle's world line (see figure 4). Note that the retarded field is generated by portion of 
the world line ( that corresponds to the particle's history before t ret {y). Since the stress-energy 
tensor is quadratic in field strengths, we should twice integrate it over £. There are also two 
variables which parametrize S t . In curvilinear coordinates (t, t\, t2, s, ip) the surface integral 
(4.1) becomes 



Pen 



[ da T 0a (4.1.4) 

t h fc? 2tt t t fc° 2tt 

J dti j dt 2 j dR j dp>Jt\% + J dti J dt 2 J dR J dtpJt\ 



Ma 
l 12 



with Jacobian 



■J - (l-g||cos^)ff (i.i.o) 
da 

+ q dR C ° S * I K 
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Figure 4: Outgoing electromagnetic waves generated by the portion of the world line that corresponds 
to the interval — oo < ti < t\ combine within the gray disk with radius k\ = t — t\. Their contribution is 
given by the first fourfold integral in eq.(4.1.4). If the domain of integration t\ < ti < t the waves joint 
together inside the dark disk with radius k% = t — t' 2 . The second fourfold integral in eq.(4.1.4) describes 
them. 

The integrand 

2v< = /#/g )A - l^ffttfJ (4-1.6) 
describes the combination of field strength densities at y £ Y, t 

generated by emission points z(t±) G C and .z(i 2 ) £ C- Symbol c a denotes the factor + (K a -v a ) 
involved in eq.(4.1.1). 

The first multiple integral calculates the interference of the disk emanated by fixed point 
z(ti) £ ( with radiation generated by all the points before the instant t\. The second fourfold 
integral gives the contribution of points after t\ (see figure 4) . 

It is worth noting that time variables t\ and parametrize the same world line (. Coordinate 
transformation (4.1.3) is invariant with respect to the following reciprocity: 

T : t\ <-> t2, a <-> (3, <p 1— > if + it. (4.1.8) 

This symmetry provides identity of domains of fourfold integrals in energy-momentum (4.1.4). 

It is obvious that the support of double integral j 1 ^ dt\ (U2 coincides with the support 

of the integral /I dt2 J*! 2 -*, dt\. Since instants t\ and £2 label different points at the same world 
line C one can interchanges the indices "first" and "second" in the second fourfold integral 
of eq.(4.1.4). Via interchanging of these indices we finally obtain J^^dti dt<i instead of 
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initial dt\ f£ dt 2 - Taking into account these circumstances in the expression (4.1.4) for 
energy-momentum carried by electromagnetic field we finally obtain 



2?r 



f cm 



t h 

dh J (It - 2 J dR J (/,.:■./ 1" 

oo — oo 



Oa 



where new stress-energy tensor is symmetric in the pair of indices 1 and 2: 
27rt 0a = 2tt (t?£ + tij?) 

_ f OA fa , f OA f a 1 Oa [ r fiv r (2) , r p.v n (1) 

- /(I)/ (2)A + / (2)/ (1)A ~ jV J(l)hv + J(2)J^ 



(4.1.9) 



(4.1.10) 



4.2. Angular integration of energy-momentum tensor density 

We see that it is sufficient to consider the situation when t\ > t 2 . The smaller disk C\{z\, t— ti) C 
C 2 (z 2 ,t — t2) is filled up by nonconcentric circles with radii R 6 [0,/cj 1 ] (see figures 4 and 5). 
To calculate the total flows (4.1.9) of electromagnetic field energy and momentum which flow 
across the plane Tit we should integrate the Maxwell energy-momentum tensor density (4.1.10) 
over angular variable <p, over radius R and, finally, over time variables t\ and t 2 . Integration 
over ip is not a trivial matter. The difficulty resides mostly with norms ||i^ a || 2 = —^apK^Ka of 
separation vectors K a = y — z a which result in elliptic integrals. To avoid dealing with them we 
modify the coordinate transformation (4.1.3). We fix the parameter (3 in such a way that the 
norm ||i^i|| 2 becomes proportional to the norm ||i^2|| 2 : 



\ Kl f = -^\\K 2 \\\ 
a 



(4.2.1) 



Keeping in mind identity a + f3 = 1, we arrive at the quadratic algebraic equation on (3 which 
does not contain the angle variable: 



R 2 = a (k^) 2 + I3{k° 2 f - a(3cf. 
We choose the root which vanishes when R = k\: 



(4.2.2) 



(-(k° 2 ) 2 + (fc?) 2 + q 2 + VB) D= [(klf - (fe?) 2 - q 2 ] 2 - 4q 2 [(A:?) 2 - R 2 } . 

(4.2.3) 

If q 2 tends to zero while t\ 7^ t 2 , it becomes the unique root of the linear equation on (3 originated 
from eq. (4.2.2) with q 2 = 0. 

If R = the R— circle reduces to point A with coordinates (z\ — Poq 1 ) where (3q = (3\ R=0 - If 
R = k® then (3 = and the circle is centered at z l {t\) (see figure 5). 

Changing the variable of integration from R to (3 transforms two inner integrals in the fourfold 
integral (4.1.9) as follows: 

fc? 2tt 2tt 

dR Jd<p(l- <?|§cos^) R = Jdf3jdip - qRcos^ . (4.2.4) 

p 
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Figure 5: The interference picture in a plane S t . The points z{t\) e C an d z (t2) € C, *2 < ti, emanate 
the radiation which filled up the disks centered at z\ and 22, respectively. The gray disk with radius 
k\ = t — ti is filled up by nonconcentric circles centered at the line crossing both the points z\ and z-i- 
If parameter vanishes the circle is centered at z\\ its radius is equal to k\. If = 0q < the circle 
reduces to the point A labeled by the box symbol. In case of intermediate value 0o < < we have the 
circle of radius R with center at point O between Z\ and A. 

Having differentiated eq. (4.2.2) with respect to we obtain new Jacobian 



It is straightforward to substitute the field densities (4.1.7) evaluated at instants t\ and ti 
into expression (4.1.10) to calculate the integrand of the multiple integral (4.1.9). It is of great 
importance that the square of norm ||K a || is proportional to J (see eqs.(B.6) derived in Appendix 
B). For the first term, Jt°^ , we obtain the following cumbersome expression 



J = (1/2) [(k° 2 ) 2 - (fc?) 2 - q 2 ] +/3q 2 - qRcostp. 



(4.2.5) 



jt: 



12 — 



e -I't'(vf4) - {/f°(A) - I'f J (A„)} . 




) 



} 



(4.2.6) 



The functions 



A = a 



d 2 a da da 



A = 



(4.2.7) 



dt\dt2 dt\ dt2 ' 



dhdt 2 
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depend on Synge's world function (3.2) of two timelike related points, z{t\) G C and zfo) G (, 
taken with opposite sign. Symbols / and I' denote /3-dependent factors 



t = + (4.2.8) 



s/=J5a' V a V -13 

Each second order differential operator 

T a = D a 7 r^- + B a -^- + C a ^- + A a (4.2.9) 

Ot\dt2 Oti Ot2 

has been labeled according to its dependence on the combination of components of the separation 
vectors K\ and K<i pictured in figure 6 (see Appendix B) or on Jacobian (4.2.5). The components 
of these vectors are involved in (/2-dependent coefficients 

D a = J--?-, B^ = ^^- 2 (4.2.10) 



27r(r 1 ) 2 r 2 ' 27r(n) 2 (r 2 ) 2 

where factor a is replaced by K?K$, Kf, K$, J or 1 for if/, ff , if, f J or f°, respectively. 

Integration of the electromagnetic field's stress-energy tensor over the angular variable is the 
key to the problem. All the (/^-dependent constructions are concentrated in the coefficients of 
differential operators of the type in eq. (4.2.9). We introduce a new operator 

dt\dt2 at i ot2 

where the script letters denote the coefficients (4.2.10) integrated over (p. (The integration is 
performed in Appendix C). 

To distinguish the partial derivatives in time variables we rewrite the operator (4.2.11) as 
the sum of the second-order differential operator 

d 2 d ( dV a \ d ( &D a \ 

r -^^ + 5h( B, --gS-) + 5s( c '-«r) <4212) 

and the tail 

. d 2 V a dB a dC a , Jni „ 

OT1OT2 ut\ Ot 2 

For a smooth function /(ti , ^2) we have 

f a (f) = fl a (f) + f7r a . (4.2.14) 
In Appendix C we derive the relations 

7T° = 0, vr J = (4.2.15) 



^7' n > =v i{ L -w) 

8V J \ a j */, 8V J 
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which allow us to rewrite the integral of eq. (4.2.6) over ip in terms of differential operators IP: 



p2tv 

Jo ' 



d(pJt^2 



fr a/3 

ii 12 



d_ 
dh 
_d_ 

dh 



d 2 



a 



dhdh 
dhdt 2 



+ m 



da \ 



dt 2 [ ° 2Vl 

d 

di 2 



8 2 a 
dhdt 2 

DVlV2 dh) 



e -^{m°(x)-i'fi J (x )}. 



(4.2.16) 



Since the operator in eq.(4.2.12) is the combination of partial derivatives in time variables, the 
angular integration gives the key to the problem. 

Setting a = and /? = i in eq.(4.2.16) we obtain the first term of the mixed space-time 
components of the stress-energy tensor (4.1.10). We add the term where indices 1 and 2 are 
interchanged. Since zeroth components of the separation three- vectors K\ and K 2 do not depend 
on ip, the final expression get simplified: 



2tt 



dipJt 0i = 



d_ 



d\ 2 
dh 



-i'fi J {v\ + v\) 



~dU 



+ n° (4Ai + ^a 2 ) 



(4.2.17) 



d_ 

dh 



where 



, n da 
dh 



, n 9a 

Oto 



(4.2.18) 



Similarly we derive the zeroth component t . Setting a = and (3 = in eq.(4.2.16) we 
obtain the first term of energy density. Since it is symmetric in the pair of indices 1 and 2, the 
second term, t%l, doubles it. The integral of energy density t 00 over the angular variable has the 
form 



dpJt 00 = e 2 \lU (K) - l'tl J (p) 
Jo L 



where 



k = 



(i 



d 2 a o da da 1 da da 

m^ 2 +kl W 1 +k2 dT 2 ~2dhdT 2 +afI 



(4.2.19) 



(4.2.20) 



1 d 2 c 



2 dhdt 2 



+ 1. 



We see that the integration of the electromagnetic field's stress-energy tensor over (p yields 
integrals that are functions of the end points only. In the next subsection we classify them and 
consider the problem of integration over the remaining variables. 
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4.3. Integration over time variables and (3 



Our purpose in this section is to develop the mathematical tools required in a surface integration 
of energy-momentum tensor density in 2+1 electrodynamics. Integration over angle variable 
results the combination of partial derivatives in time variables: 

ft rti 



P£n(0= "7T 

2 



/ dh [ 

J —OO J — 

f dt 2 f 

J -co Jti 



dt 2 



d 2 Gf 2 dGf dG^ 



dt\dt 2 dt\ 



By virtue of the equality 

r° dG(p,t u t2) 

ho dt a 



d_ 



00 



d(3G((3,h,t 2 ) 



+ G(p ,h,t 2 ) 



du 



dp (h,h) 
dt n 



(4.3.1) 



(4.3.2) 



the triple integral (4.3.1) can be rewritten as follows: 



Pcm(t) 



lim 



o 



d(3G 



— e 



+ e- 



/. 



k°->0Ji3 
t 

dt 2 lim 

t 



12 



*2=t 

— CO 

dG% 



dt 2 lim 

oo fe?^0 _ 



(4.3.3) 



d/3 



Ho 



du 



+ e" J dt 2 lim 

k°=k°-At J-co fc?^0 



d/3G? 



dt\ lim 

At^O 



+ 



/ dh [ 

J —CO J — 

1 f dt 2 [ 

J —oo JU 



dt 2 



t 



k°=k°+At 



j —( 



dt\ lim 

t 2 ^-co 



[°dpG% 
UPo 



-dG\ 2 



dti 



dt 



J P=Po 



d ( 3 ra \ 9 ^0 



In Appendix C we calculate the functions 
G° 12 = IV°k-I'V j h 



(4.3.4) 



G? 



-/3 v 

K 



a 



P 



dp VIM 



a . 2 w - 



-/3'lr 2 113 



a 



9/3 V lira 



G 



12 



G\ 



G% 



I 

2 

2 



dh 



(4.3.5) 



9Ai 



dA 



^ (aq^j + rM) + ^ (-PqWj + r° V j) + (t;*A 2 + ^i) v? 



9 



2^ + ^) a/3 



2 |H| 3 



9^ 



[aq^ 2 2 + r 2 %*) + ^ (-/3<f v| + + {v\\ 2 + v * Ai) 



1 + «£) 



9 



0/3 V ||r 2 
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involved in these integrals. 

1°. Integrals where ti — ► t. Equality (B.10) implies that the lower limit (3q tends 
to if = t — t\ vanishes. The upper limit is equal to zero too. Then the integral over 
parameter [3 vanishes whenever an expression under integral sign is smooth. So, we must limit 
our computations to the singular terms only. They are performed in Appendix D; these integrals 
do not contribute in the energy-momentum at all. 

2°. Integrals where ti = ti- According to eq.(B.lO), the equality t\ = ti yields sin$o = 1 
and lower limit (3q = — tan^o tends to — oo. The small parameter is the positively valued 
difference At = t\ — ti- The integration is performed in Appendix E; the resulting terms belong 
to the bound part of energy-momentum (to that which is permanently "attached" to the charge 
and is carried along with it.) 

3°. Integrals where t2 — > — oo. Equality (B.10) implies that the lower limit (3q tends to 
if = t — ti increases extremely. Then the integral over parameter [3 vanishes whenever an 
expression under integral sign is smooth. So, we must limit our computations to the singular 
terms only. They are performed in Appendix E; the resulting terms belong to the bound 
electromagnetic "cloud" which can not be separated from the charged particle. 

4°. Integrals at point where (3 = (3q. In this case the radius of the smallest circle 
pictured in figure 5 vanishes and it reduces to the point A. The contribution in p^ m is given 
by the last line of eq. (4.3.3). In Appendix F we present the integrand as the combination of 
partial derivatives in time variables and nonderivative tail. After integration over t\ or ti, the 
derivatives are coupled with bound terms obtained in Appendix E; the sum is absorbed by 
three-momentum of bare particle within renormalization procedure. The tail contains radiative 
terms which detach themselves from the charge and lead independent existence. 

Summing up all the contributions 2° — 4° we finally obtain 



The finite terms which depend on the end points only are noncovariant. They express the 
"deformation" of electromagnetic cloud due to the choice of coordinate-dependent hole around 
the particle in the integration surface We neglect these structureless terms. The single 
integrals describe the covariant singular part of energy-momentum carried by electromagnetic 
field. The first term in between the square brackets of eqs. (4.3.6) and (4.3.7) cannot be rewritten 
as the partial derivative in t\ or ti- It determines the radiation reaction in 2 + 1 electrodynamics. 




(4.3.6) 




ti=t 



(4.3.7) 
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5 Conclusions 



In the present paper, we calculate the total flows of (retarded) electromagnetic field energy, 
momentum and angular momentum which flow across the plane £$ = {y G M3 : y° = t}. The 
field is generated by a point-like electric charge arbitrarily moving in flat space-time of three 
dimensions. The computation is not a trivial matter, since the Maxwell energy-momentum 
tensor density evaluated at field point y G £i is nonlocal. In odd dimensions the retarded field 
is generated by the portion of the world line £ that corresponds to the particle's history before 
t ret (y). Since the stress-energy tensor is quadratic in field strengths, we should twice integrate it 
over £. We integrate it also over two variables which parametrize T, t in order to calculate energy- 
momentum and angular momentum which flow across this plane. Thanks to the integration we 
reduce the support of the retarded and advanced Green's functions to particle's trajectory. 

The Dirac scheme which manipulates fields on the world line only is the key point of inves- 
tigation. By fields we mean the convolution of three-velocity and nonlocal part of the retarded 
strength tensor evaluated at point z(t\) G C; the torque of this "Lorentz 6>-force" arises in elec- 
tromagnetic field's total angular momentum. (The singular 5-term (2.7) is defined on the light 
cone; it is meaningless since both the field point, z(t±), and the emission point, zfe), lie on 
the time-like world line). The retarded and the advanced quantities arise naturally. The re- 
tarded Lorentz self-force, as well as its torque, contains integration over the portion of the world 
line that corresponds to the interval —00 < t<i < t\. Domain of integration of their advanced 
counterparts corresponds to the interval t\ < £2 < t. 

Noether quantity G^ m carried by the electromagnetic field consists of terms of two quite 
different types: (i) singular, Gg, which is permanently attached to the source and carried along 
with it, and (ii) radiative, G^, which detaches itself from the charge and leads independent 
existence. The former is the half-sum of retarded and advanced expressions, integrated over £, 
while the latter is the integral of one-half of the retarded quantity minus one-half of the advanced 
one. Within regularization procedure the bound terms Gg are coupled with energy-momentum 
and angular momentum of the bare source, so that already renormalized characteristics Gp art 
of charged particle are proclaimed to be finite. Noether quantities which are properly conserved 
become 

U — ^part + ^R- 

The regularization procedure which relies on energy-momentum and angular momentum bal- 
ance equations is proposed. Energy-momentum balance equations define the change of particle's 
three-momentum under the influence of an external electromagnetic field where loss of energy due 
to radiation is taken into account. The angular momentum balance equations explain how this 
already renormalized three-momentum depends on particle's individual characteristics. They 
constitute the system of three linear equations in three components of particle's momentum. Its 
rank is equal to 2, so that arbitrary scalar function arises naturally. It can be interpreted as a 
dynamical mass of dressed charge which is proclaimed to be finite. A surprising feature is that 
this mass depends on the particle's history before the instant of observation when the charge is 
accelerated. Already renormalized particle's momentum contains, apart from the usual velocity 
term, also nonlocal contribution from point-like particle's electromagnetic field. 

Having substituted this expression in the energy-momentum balance equations we derive a 
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three-dimensional analogue of the Lorentz-Dirac equation 



e 2 [ T ds 

ma» = eUr, a F£ (r) + -aM — + eu T , a F^ t . 

1 J-oo ^2a(T, s) 

The loss of energy due to radiation is determined by work done by Lorentz force of point-like 
charge acting upon itself. The nonlocal term which is proportional to the particle's acceleration 
provides finiteness of the self-action. The third term describes the influence of an external field. 

In this paper we develop a convenient technique which allows us to integrate nonlocal stress- 
energy tensor over the spacelike plane. The next step will be to implement this technique to a 
point particle coupled to massive scalar field following an arbitrary trajectory on a flat space- 
time. The Klein-Gordon field generated by the scalar charge holds energy near the particle. 
This circumstance makes the procedure of decomposition of the energy-momentum into bound 
and radiative parts unclear. 

In Ref. [20] the remarkable correspondence is established between dynamical equations 
which govern behaviour of superfluid 4 He films and Maxwell equations for electrodynamics in 
2+1 dimensions (see also refs.[21, 22]) 2 . Perhaps the effective equation of motion (3.30) will be 
useful in study of phenomena in superfluid dynamics which correspond to the radiation friction 
in 2 + 1 electrodynamics. 
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Appendix A Uniformly moving charge (conserved quantities) 

The simplest field is generated by an unmoved charge placed at the coordinate origin. Setting 
z = (t, 0, 0) and u = (1, 0, 0) in eq.(2.11), one can derive that the only nontrivial components of 
static field are: 

y°-r 

^ f dt y l . . „ . 

F " ~ * J TCT^f ^y ^ 



= —e- 



V(y°-t) 2 -r 2 



r 2 y u — t 



t=y°-r 

= e— 



r := \J (y 1 ) 2 + (y 2 ) 2 is the distance to the charge. Having performed Poincare transformation, 
the combination of translation and Lorentz transformation, we find the field generated by a 
uniformly moving charge: 

Uak/3 - Upka 

F a /3 = e — . (A.2) 



2 I wish to thank O.Derzhko for drawing these papers to my attention. 
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Here r = —T] a p(y a — z a {s))u^ is the retarded distance where the particle's position z a (r) = 
Zq + u a r is referred to the retarded instant of time s. We denote k a the null vector y a — z a {s) 
rescaled by the retarded distance, i.e. 

k- = "°-^">. (A.3) 
r 

It is straightforward to substitute eq.(A.2) in eq.(4.2) to calculate the electromagnetic field's 
stress-energy tensor: 

2vrT a/3 = ^ (u a kP + vPk a - k a kP + \*f^ . (A.4) 

Now we calculate the electromagnetic field momentum 

Pem(t) = [ da,T^ (A.5) 

where an integration surface S t is a surface of constant y°. We start with coordinate transfor- 
mation (y° ,y l -,y 2 ) ► (r,s,<p) locally given by 

y a = z a {s) + rk a , k a = A a a ,n a ' (A.6) 

where n a = (1, cos (p, sin y?). The Lorentz matrix A determines the transformation to the parti- 
cle's comoving Lorentz frame where the particle is at rest. To adopt these curvilinear coordinates 
to the integration surface H t = {y G M3 : y° = t} we replace the retarded distance r by the 
expression 

r = *-w (A - 7) 

where t is the observation time. On rearrangement, the final coordinate transformation 
(y°, y 1 , y 2 ) 1— > (t, s, ip) looks as follows: 

y° = t, y W(s) + (t-s)^. (A.8) 



Differentiation of this coordinate transformation yields the differential chart 

1 



g t = ^{k% + k%) (A.9) 



k 



e v - (* s ) y fc o (kP) 2 ^ j e ' 1 

where k^ = dk a /dip and v i = 7 _1 u\ Their scalar products are the components of metric tensor 
g of Minkowski space M3 as it is viewed in curvilinear coordinates (A.8). To calculate the 
determinant of g it is sufficient to know that 

gtt = 0, gtip = 0, g ts = — p-, g w = ^ 2 . (A.10) 
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The surface element is given by da$ = ^/—gdsdip where 

^? = 7-^| (A.11) 

is the Jacobian of coordinate transformation (A. 8). Electromagnetic field momentum (A. 5) takes 
the form: 

t 2w 

p? m = ^ j dsj d<p^- (uV + uW - k°k? + (1/2)^) . (A.12) 

-oo 

The angular integration can be handled via the relations 

J dipk a = 2jru a , J dpk a k? = toru a vP + 7T7/ Q/3 . (A.13) 



After trivial calculations we arrive at the logarithmic divergence 

= (A - 14) 

— oo 

as could be expected for the two-dimensional Coulomb potential. 
Appendix B Coordinate system 

The coordinate transformation (4.1.3) is associated with two points, z(t\) and zfa), on an 
accelerated world line ( (see figure 6). Its differentiation yields differential chart 

et = e - Q-^ e i' 
e R = n J e f - 

e v = Rn^eji (B.l) 
where n = (cos tp, sin p), = (— sin ip, cos (p) and ey = e^y. The orthogonal matrix 

/ nl —n 2 n \ 

v=[ * / (B.2) 

where n* = q % jq rotates space axes until a new x— axis be directed along two- vector q := 

z(ti) — z(i2) (we denote q := \A?)- I n new coordinates, three-vectors K a = y — z(t a ),a = 1,2, 
pointing from points of emanation z{t a ) := z a to an observation point y G S t (see figure 6) have 
the following forms: 

K° a = t-t a :=k° a , K = Jik{ (B.3) 

where 

k\ = -(3q + Rcosp, kf = Rsin<p (B.4) 
k\ = aq + Rcosip, k\ = Rsinip. 



25 




Figure 6: The separation vector K a is a vector pointing from point of emission z(t a ) := z a to point of 
observation P e S t with coordinates (y°, y 1 , y 2 ). The integrand t 0a involves also one-half of square of 
three-vector q := z\ — zi = if 2 — K\ (double Synge's function) and its partial derivatives with respect to 
time variables. Space components of q determine the orthogonal matrix (B.2). 

The norms ||-ff a || = \/—{K a ■ K a ) of the separation vectors K\ and K2 pictured in figure 6 
contain the angular variable: 

-{Kx-Kx) = (fc?) 2 -/3 2 q 2 + 2(3qRcos(p - R 2 , 

-(K 2 -K 2 ) = {klf -a 2 q 2 - 2aqRcosip - R 2 . (B.5) 

Substituting the right-hand side of eq. (4.2.2) for R 2 in these expressions and comparing them 
with Jacobian (4.2.5) leads to the important relations: 

||#i|| 2 = -2/3J, \\K 2 \\ 2 = 2aJ. (B.6) 



They immediately follow: 

J 1 



(B.7) 



To concretely compute integrals over (3 we clarify the mathematical sense of this parameter. 
Since eq.(B.6), we parametrize the ratio ||ifi||/||if2|| by angle variable #: 



si ni ? = A/^. (B.8) 
V a 

Having inserted it in expression (4.2.2) we obtain algebraic equation on sin 2 1?. Its solution is as 
follows: 
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where D is defined by eq. (4.2.3). If R = then i9 = and parameter [3 vanishes. If R = we 
obtain the lower limit of the integral over (3. We denote it as (3q = — tan 2 #o where 



sintf = (v^S- V2a) (B.10) 

and 

S = 5 (*§ + fe?) 2 - iq 2 , a = 1 (A; 2 ° - fe?) 2 - Iq 2 . (B.ll) 



Appendix C Integration over angular variable 



To calculate the energy- momentum (4.1.9) carried by the electromagnetic field we should perform 
the integration over angle first. When facing this problem it is convenient to mark out ip- 
dependent terms in expressions under the integral sign. In the Maxwell energy-momentum tensor 
density we distinguish the second-order differential operator T a with ^-dependent coefficients 
(see eq. (4.2.9)). Having integrated this operator over tp we obtain the operator l~ a which can be 
decomposed into a combination of partial derivatives in time variables Il a given by eq.(4.2.12) 
and tail 7r a of the type in eq.(4.2.13). 

This Appendix is concerned with the computation of the tails. Equipped with them we 
express the integrand as a combination of partial derivatives in t\ and t 2 . It allows us to 
integrate the electromagnetic field's stress-energy tensor over the time variables, as well as over 
P. 

To implement this strategy we must first integrate the coefficients (4.2.10) over the angle 
variable. We start with the simplest one 

2vr J ~ r nr 2 



where (/^-dependent numerator a is equal to the Jacobian (4.2.5) or to 1. Our task is to rewrite 
the integrand as a sum of term with denominator n and term with denominator r 2 . To do it 
we introduce a new layer of mathematical formalism and develop convenient technique. 

We introduce null vector n = (1, cos cp, sirup) which belongs to the vector space, say V, 
such that i , ii, and i 2 is its linear basis. We shall use rj a p = diag(— 1,1,1) and its inverse 
rf-P = diag(— 1, 1, 1) to lower and raise indices, respectively. We introduce the pairing 

(•) : VxV^R (C.2) 

(a • b) » Va pa a bP (C.3) 

which will be called the "scalar product" . 
We express the ^-dependent constructions 

r a = -(K a -v a ), c a = 7 ~ 2 + (K a ■ v a ) (C.4) 

as the scalar products — (r a • n) and (c a • n), respectively. We shall use sans-serif letters for the 
components of timelike three-vectors r a G V and c a € V 

r° = fc?+^(viq), r° = fc 2 ° - a(v 2 q), r{ = Rv a W; (C.5) 
c? = -7f 2 + /3(viq), = -7 2 ~ 2 -a(v 2 q), d> a = Rvjui j . 
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The Jacobian (4.2.5) becomes the scalar product (J • n) = J + Ji cos 92 + J 2 sin tp, where 

Jo = /3q 2 + (1/2) [(k° 2 ) 2 - (fc?) 2 - q 2 ] , h = -qR, h = 0. (C.6) 

We introduce the dual space of one- forms, say W, with basis u>°, cj 1 and a) 2 such that a) M (i^) = 
by where io, ii, \i constitute the basis of V. The wedge product L = a A b of two one forms a 
and b constitutes two-form 

L = (ao&i — ai&o) cj° A ci 1 + (ao&2 — «2^o) w° A w 2 + («i&2 — a2&i) cj 1 A u> 2 . (C.7) 

We introduce dual three-vector L =*L with components 

L« = (C8) 

e a /37 denotes the Ricci symbol in three dimensions: 

{1 when a/37 is an even permutation of 0,1,2 
— 1 when a/37 i s an °dd permutation of 0,1,2 (C9) 
otherwise. 

We raise indices in eq.(C.8) and define the vector product of two vectors, a and b: 

l a = e a ^b v . (CIO) 

Tensor 

e\„ = e a ^ Wlu (C.ll) 

has the components 

e\ v = e ^, e\ v = -e 1 ^, e\ v = -e 2 ^ . (C.12) 
It is interesting that tensor 

e\^u = e^VaWPuVju (C.13) 
is equal to e XfxI/ taken with opposite sign. 

Now we calculate the double vector product 

[A[BC]] a = e a ^e\ v WC. (C.14) 

Since 

£ Q /3 7 £V = -^vftfSv + tPvVPii (C.15) 

we arrive to the unusual rule 

[A[BC]] = -B(A-C) + C(A-B) (C.16) 
instead of the well-known law acting in space with Euclidean metric. 
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To simplify the denominator r±r 2 in the integrand of eq.(C.l) as much as possible we rewrite 
27r-periodic functions (C.4) as follows: 



r a = -r a ,o- Pa sm((p + 4> a ), p a = ^/r^ + r 2 2 . (C.17) 

(We recall that r a is the scalar products (r a • n) taken with opposite sign, components are 
given by eqs.(C.5).) Shift in argument of a-th function is determined by the relations 

r a ,i = Pa sin 4> a , r aj2 = p a cos<f> a . (C.18) 

After some algebra we rewrite the integrand of eq.(C.l) as the following sum 

a _A\ 2 - CVi cos((^ + fa) + A^ + Cf 2 p 2 cos(y + 2 ) 1 



ri r 2 ri r2 

Coefficients Af 2 , A 21 and Cf 2 satisfy the vector equation 

-A a 12 r 2 -A a 21 r 1 + C? 2 L 12 = a (C.20) 

where by ri and r 2 we mean three- vectors with components in eq.(C.5) and L12 = [fi^]- 

To solve equation (C.20) we postmultiply it on the vector product [ril_i2], then on the vector 
product [r2l_2i] and, finally, on L12. After some algebra we obtain 

,q _ (t ar l] • L 12) 4 q _ M • Ul) r a _ (g ; j-lg) , p „ n 

-^12 -^21 -^12 

where the denominator D\ 2 = (L\ 2 ■ L\ 2 ) is symmetric in its indices. 
Substituting eq.(C.19) into eq.(C.l) and using the identities 



1 f 2w dp 



27r Jo r a - p a sm{(p + (j) a ) 
-' lT ' cos(99 + a ) 



(C.22) 



1 f Z7r , cos( 

— / dip-r; 

2vr Jo - p a 1 



,sin(ip + (p a ) 

yields 



l r l|| ll r 2| 

after integration over <p. 

Now we turn to the calculation of the coefficient 

1 f 2n ac 2 



A a A a 



B a = — dip-f^-. (C.24) 
2itJ ri(r 2 ) 2 V ; 

Equipped with the relations in eq.(C.19) we rewrite the integrand as a sum of terms which are 
proportional to the 1/ri, l/r 2 , and l/(r2) 2 , respectively. Using the identities 



dtp r;; 



1 /" Z7r 



2vr7 [r°-p o sin(p + a )] 2 

1 C0S((^ + (j) a 



(C.25) 



d(p rr = 

2W [rO- Pa sin(^ + ^)] 2 
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and taking into account the relations in eq.(C.22) gives 



r 2 



1 (a ■ r 2 )(c 2 • r 2 ) 1 
-\J2 ■ rij + 



+ 



9 4 a A C2 — 



£>2i || 
([ari] • [c 2 ri])" 



4 C 2 4 a _l A a 4 C2 



12^21 



(a • c 2 )(ri • r 2 ) 



£>12 



The resulting expression for the term 



1 f 2n 
2Wo 



aci 



(C.26) 



(C.27) 



(n) 2 r 2 

can be obtained by interchanging indices 1 and 2 in the right-hand side of eq.(C.26). 

After a routine computation based on the repeated usage of relation (C.19) we find the last 



term 



A a — — 



2?r 



dip 



ac x c 2 

(n) 2 (r 2 y 



(C.28) 



B12 , B 21 (a ■ ri ) 
Ti — H" + Ti — iT + — ' l — Ti — 11T 

M ll r 2ll Ki 



(a • r 2 ) 



l r 2| 



where 



7 _ o^ci^ca ([ciri] • [c 2 ri]) 



(C.29) 



B 



12 



o Aa aci A c 2 1 941 4 C 2 4C1 1 94C1 /IC2 4a 1 Aa 
— ° /1 12^ 1 12^ 1 21 ^ 0/l 1 2 yl 1 9^91 ' 19 19 91 T^ 1 !' 



12^12^21 



l 12^ l 12^ l 21 



l 12 



([ciri] • [c 2 r 2 ]) ([cir 2 ] • [c 2 ri]) 



A. 
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£>12 

C1 ([c 2 ri] • [ari]) ([ciTi] • [ari]) ([c 2 ri] • [ar 2 ]) ([ciri] • [ar 2 ]) 

^12 ^12 ^12 -^12 



+ 
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and the others, -B21 and J 2 , can be obtained by interchanging indices 1 and 2. 

We now turn to the differentiation of coefficient (C.23) with respect to time variables t\ 
and t 2 . We will use Latin indices a and b which run from 1 to 2 (a 7^ b). We introduce new 
denotations k\ = a and k 2 = (5 for time-independent variables (3 and a = 1 — (3. Differentiation 
of V a is based on the relations obtained via differentiation of zeroth components (C.5), (C.6) 
and square of radius (4.2.2): 



dt a 
dt 



a*ai 



— = -(-l)°rS-«a(v a q), 



-2 Ka r°. 



They immediately give 
d(r a • r a ) 

<9(r fe • r a ) 



= 2(r a • c a ), 



<9(r a • r a ) 2K b 



5t b 



dt b R 2 
(r a • Cft) + [r°(r b • r 6 ) - r£(r a • r b )] 



[r°(r a • r b ) - r°(r a • r B )] 



(C.30) 
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and eventually gives 

d ( 1 \ _ (Co ■ r a ) d ( 1 \ n b r°(r a ■ r b ) - r°(r a • r a ) 



dta \\K\\J \K\\ 3 ' dt b \ \\r a \\J ||r a || 3 R 2 
After some algebra we also obtain 



(C.31) 



= (-i r (r a .r a ) + (J.c a ) + ^[j (r a .r a )-r°(J.r a )] (C.32) 
= (-l) b (r 6 -O + §[J°(raT b ) + r°(J.r fe )-2r0(j.r a )]. 
Usage of these relations allows us to calculate the following derivatives 

^ = ^ + ^X+ ([Jr ^ [Car6]) (C33) 



9t a a0 fta aft oa D ab 

- b + 2A c * h A J ab - ([Jra] p [C&ra]) " § (J° + M 2 + MO 



where latin indices a and 6 run from 1 to 2, a ^ b. 

Having differentiated eq.(C.23), after a straightforward calculations we derive the following 



relations 

^- = B J -i-(-^-r). (C.34) 



d pJ „j d f a \ dV J j d ( p 



dti 0/3 V II r 2 1| y ' dt 2 d(3 

Further we find out the expression dC J /dt 2 and prove the identity 



9t 2 9ti V ^ 

(One can derive dB J /dt\, subtract it from „4 J and compare the result with d/dt 2 {C J -&D J /dh).) 

Similarly we derive an analogous equality where Jacobian J with components (C.6) is replaced 
by unit three-vector o = (—1,0,0). Having substituted o for a in the expressions (C.23), (C.26) 
and (C.28) we obtain the terms T>°, B° and -4°, respectively. The remaining term, C , can be 
obtained from £>° via reciprocity. The derivatives of coefficients A° ab are as follows: 



9 Kb ,0 A c a | ,c ,0 | ([° r a] • [Cqr&]) 



fli " ab-^ba ' ab ba ' t-j 



ab 



(C.36) 



9A ab _ nACb A ([or a ] • [CfcTq]) Kf, , .0^ 

- 2^ ab A ab — ^(o + r 2 A 12 + h A 21 ). 



Substituting these into equality 



dv° d ( A\ 2 , ^ 



9t a 9t a VII rill || r 2 1 
and using the identities (C.31) yields 



+ tt^t: (C.37) 
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And, finally, we calculate the partial derivative dC°/dt 2 , subtract it from A and compare the 
result with d/dh(B° - &D°/dt 2 ). We obtain 



Now, we calculate the tail 



where 



^ =A *-^-m; + m^r 2 (c - 40) 



K = 7T / d^- 2 -, = «- / (C41) 

r2n r ,2k K*C X C 2 



1 [ , : - K 2a 4 «_j_r 

- - 2vr7o Vi) 2 r 2 ' " " 2vr i 



(ri) 2 (r 2 ) 2 ' 



The zeroth component, = does not depend on ip. Inserting relations = k^D°, 
Bl = k°B°, C° = /c°C°, and A° a = k°A° into eq.(C40) and taking into account identity (C.39) 
yields 

-flfe' 7r2 = c -a*r- (c - 42) 

Space components, if*, depend on ip. They can be expressed as the scalar product (K l a ■ n) 
where components of three-vectors K l a G V are as follows: 

K* 1)0 = -M K 2j0 = a^, K^iM, Kl >2 = W 2 , (C.43) 

where w*j are components of the orthogonal matrix (B.2). Having substituted K* for a in 
expressions (C.23), (C.26), and (C.28) we obtain the terms V l a , B l a and A l a , respectively. The 
last term, C l a , can be obtained from B l a via reciprocity. To differentiate them we need the 
equalities 

= (K ■ ca) - vtf - § [^(ra ■ r B ) + r°(K* a ■ r a )] 
= -^ + g[r°(Ki.r a )-Ki !0 (r 6 .r a )-2r^(Ki.r 6 )] 

= (Ki-c b )-g[Ki )0 (r b .r b ) + r 6 °(Ki.r 6 )] 
3(K* -r a ) _ K fe r n 



^ = ^[ r a(KL-r fe )-K* i0 (r 6 .r a )-2r°(K* -r a ) 



in addition to eqs.(C.30) and (C.31). 
The derivation of equalities 



= vlV° + ^(-(3 q Wl + r ° 2 vl) (C.44) 
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dt 2 ||ri| 



= W Kvl + r^) 



is virtually identical to that presented above, and we shall not bother with the details. Finally, 
after a straightforward (but fairly lengthy) calculations we derive the following relations 

which generalize eqs.(C.42). 

In analogous way one can derive the equalities 

< c ' 46 > 

which arise in (^—integration of angular momentum carried by the electromagnetic field. 
We will need also the tail 

^ = S^-^-^+< (C47) 



where 



dt\dt2 dt\ dt2 



2vr J nr 2 2vr J ri(r 2 ) 2 



Cl2 " 2vri ^(n)^' " 4l2 "2 7 r7 ^ (n) W 

It can be obtained by means of covariant generalization of previous relations. Setting a = and 
(3 = and taking into account eq.(C.39) we obtain 



ati at 2 

where relations V® = k®V°, B® = k®B°, and C® = k®C° are taken into account. If a = i and 
(3 = we have 



33 



If a = and (3 = j we arrive at 



<)l. ) ( 7/ 2 
• ^ (C51) 



0*i J z dt 2 
An obvious generalization of expressions (C.49)-(C.51) is 



Appendix D Calculation of integrals where ti 



In this Appendix we compute the integrals in eq. (4.3.3) where time parameter t\ tends to the 
observation time t: 



lim / d(3G^ 2 
fc?-0 7/9b 



t 2 =t 



hm 

oo fcj— >0 



(D.l) 



dt2 lim 

oo fej— »0 L^/3q 



J0n 



Equality (B.10) implies that the lower limit f3o tends to if k® = t — 1\ vanishes. With a degree 
of accuracy sufficient for our purpose, 



(*?) S 



(k° 2 y-cf(t,t 2 y 



(D.2) 



Integrals over parameter j3 vanish whenever an expression under integral sign is smooth. So, we 
must limit our computations to the singular terms only. We expand expressions in powers of 
the small parameter (5 and then integrate them. 

Let us consider contribution p° to the energy p® m carried by electromagnetic field. Setting 

= in two-point functions (4.2.20) we obtain 



n(t,t 2 ) = v(t,t 2 )-^- + afi(t,t 2 ) 
ot 2 



(D.3) 



MM2) = - 



dt 2 



where 



it,t 2 ) = -[k° 2 + (^t)]. 



Integration of the singular part of G\ 2 given by eq. (4.3.4) gives the regular expression: 



I- /V; 2 =-Am 1+ ! vt ! r (M2) 



(D.4) 



(D.5) 
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where v t denotes the particle's velocity referred to the time of observation. According to 
eq.(D.l), we should take the function in eq.(D.5) at the end points, i.e. its value at the re- 
mote past t 2 — ► — oo should be subtracted from its value near the observation instant t. 

Similarly we calculate the integral 



lim / d(3G\ 



1 , l + |v t 
■ m 



1 + 



v t |"l-|v t |; [2a{t,t 2 )? /2 



(D.6) 



In 



1 + N 



2|v t | 1 - N/ y/2a(t,t 2 )' 
According to eq.(D.l), the result should be added to the limit 

-2K(t,t 2 ) , n(t,t 2 ) 



lim 



Gl2 dt2 



+ 



(D.7) 



and the sum should be integrated over t 2 . After a straightforward calculations we obtain the 
expression (D.5), taken with opposite sign. Their sum vanishes. Therefore pi does not contribute 
in the energy carried by electromagnetic field. 

In an analogous way we calculate contribution p\ to the momentum of the electromagnetic 
field. Integration of the singular part of function (4.3.5) over (3 gives 



lim [° dpG\ 2 



t 2 =t 



t 2 ^ — oo 



1 



2 v t 



In 



i + N 
i - N 



+ v\k\ 



VMM2) 



t 2 =t 



t-2- 



(D.8) 



It is the first term in eq.(D.l). The limit under the integral sign (second term) is as follows: 



lim 

ti-»t 



G 



12 



p=p dt 2 



We add it to the integral 



lim 

ti-»t 



d(3G\ = 1 



A) 



1 



In 



[2a(t,t 2 )f 2 df2' 
1 + N \ q l + «t& 2 <9(7 



2|v t | _ 1 

1 + N 



In 



|v t |; [2a(t,t 2 )]3/ 2 5t 2 



(D.9) 



(D.10) 



2|v t | 1 - |v t | y/2a(t,t 2 ) 
and integrate over i 2 . We arrive at the function of the end points only which annuls eq.(D.8). 



Appendix E Calculation of integrals where t 2 — > ii and t 2 — ► — 00 

In this Appendix we compute the integrals in eq. (4.3.3) where time parameters t\ and i 2 are 
equal to each other. We add also the integral evaluated at the remote past: 

(E.l) 



p%(t) = -e 2 f dt 2 lim 

7_ oc At-0 



J Pa 



dU 



+ e 



] f dti lim 

J -00 A *^o 



[°d0G% 
Upo J 



— e 



t 2 =ti-At 



t!=t 2 +At 

2 f dh lim f dpG%. 

/-OO *2^-OC Jo 
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For fixed instant t\ we assume that the limit 



lim — 



(h,t 2 ) 



t 2 



(E.2) 



is finite if the motion is infinite. (In the specific case of finite motion, A l =0.) If k® — > +00 the 
lower limit of /3-integrals 



Po 



0\ 2 



(E.3) 



1 - A 2 ' 

tends to the upper limit (zero) and we must limit our computations to the singular terms only. 

According to eq.(B.lO), the equality t\ = t 2 yields sin??o = 1 and the lower limit, Po = 
— tani?0) tends to —00. In this case the small parameter is the difference At = t± — t 2 . If the 
instants t\ and t 2 are close to each other, function /?o(t,ti,t 2 ) raises as (At) -1 : the product 
PoAt possesses finite limit. We expand the expressions under integral sign in powers of At and 
thereafter we integrate the series. The integration over (3 can be handled via the relations 



dp 



J/3 V-fia 

r° dp 
h 



ho V-P® 

f° dp 
J/3 V-Pa 



P 



= 21n(7^A) + ^) 

In (V~A) + y/a^j 



(E.4) 



-Po 



P 2 



-^Poao 



-Po 



Next we take the limit At 



«o 4 V «o 
0. Suffice it to know that 

lim At [° -4^=P 



h - m 

4 



-Po + ^fao 



lim 

At 



A, y/~Pot 
(At) 2 /° ^ 

/ft, 



ho yJ-Pa 



P Z = 



(fcO) 2 



(E.5) 



1 + VI 



We begin with the zeroth component. Setting eqs. (4.3.4) in eq.(E.l), after some algebra we 
arrive at 



(t) = -e 2 f dt 2 lim [ / *£= (^V° + K B°) 

J -00 At^O Ufa \dt 2 J 



(E.6) 



'Po 



yJ—Pot \dt 2 



+ e 



: f dh 

J —00 





Po 



ti=t 2 +At 
'0 



- / dp * I — 577— — m rdpi'4-j a 



-0H 



9/3 V ||r 2 



-, t 2 =t 1 -At 



Recall that 



1 /" 27r 
V a = — , 



rir 2 



1 /" 27r 
B a = — , 
2W0 



ac 2 



n(r 2 ) 2 



(E.7) 



where a = 1 for T> ,B and a is equal to Jacobian (4.2.5) for coefficients labeled by J. 
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The quantity r 2 = — (r 2 ■ n) can be related to the quantity r\ = — (ri • n) by Taylor expansion 
in powers of At. With a degree of accuracy sufficient for our purposes we obtain 



T2 = Ti + Atc 2 

where c 2 = (c 2 • n). Integration of (E.7) over the angular variable gives 
(a • r 2 ) 
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Expanding the integrands in eq.(E.6) in powers of At and using the relations in eq.(E.5) we 
finally obtain 
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after integration by parts of the last term in eq.(E.6). 

The first and the second terms are singular. To deal with divergences it is efficient to 
introduce the hyperbolic angles ^ and ib: 



cosh •!> = 



1*0 i 1.0 
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sinh •]/ = 
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(E.ll) 



k° - k° 

cosh ib = — - , sinh ib = 

q q 

(Function S(t, t l5 t 2 ) is introduced in Appendix B.) In these notations 
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— - [cosh(^ — ib) — l]=— sinh 2 — - - 
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[cosh(* — ib) + l] = cosh 2 — -. 



so that 
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Since the factor before the sign of limit is the total time derivative, the logarithmic divergence 
in eq.(E.lO) can be integrated by parts: 
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Taking into account that at the end points hyperbolic angles vanish, we finally obtain 
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Now we rewrite the second divergent term involved in eq.(E.lO). We expand the integrand 
in powers of At. Passing to the limit At — > 0, we arrive at: 
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Inserting these expressions in eq.(E.lO) we finally obtain 
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The calculation of momentum corrections is virtually identical to what is presented here, 
and we shall not worry with the details. It suffices to present the resulting expression 
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In the next Appendix the contribution from integrals at point (3 = (3o is found. Its bound part 
contains the terms which annihilate the divergent terms in expressions (E.17) and (E.18). 



Appendix F Calculation of integrals at point where (5 = /3q 



We now would like to extract the partial derivatives with respect to time variables from the 
integrand of the following double integral 
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Note that 
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and functions Gf 2 ,Gi, and G 2 are given by eqs. (4.3.4) and (4.3.5). 

If (3 = Po the radius R of the smallest circle pictured at figure 5 vanishes and it reduces to 
point A. Norms ||r || and \\c a \\ become zeroth component r° and c°, respectively. Hence the 
coefficients (4.2.10) get simplified, e.g. 
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In terms of two-point functions a = —l/2(q ■ q) and X = o - + /cj/c^ the angle-free functions k, \i 
and A a involved in p^ m are as follows: 
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First we set a = 0. Routine scrupulous calculations allow us to rewrite the expression under 
the integral signs in eq.(F.l) as follows: 
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This expression contains the term which is proportional to the mixed second-order partial deriva- 
tive of a which can not be rewritten as a derivative with respect to t\ or t 2 . 

In analogous way we rewrite the Pq part of space components p l em of the momentum carried 
by the electromagnetic field: 
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Having integrated expressions (F.5) and (F.6) over time variables t\ and t 2 according to the 
rule (F.l) we obtain 
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So, besides the double integral which describes the self-action which depends not only on the 
current state of motion of the particle but also on its past history, we have the integrals of 
functions of the end points only 

According to eqs.(B.ll), function £(t, t±,t2)\ t =t is equal to the function a(t, t 2 ). This cir- 
cumstance simplifies evaluation of the terms referred to this end point. We expand the terms 
near t 2 = h in powers of At = h — t 2 and take the limit At — > 0. We use the assumption in 
eq.(E.2) when t 2 — > —00. After some algebra we finally obtain: 
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Divergent terms annul their counterparts from eqs.(E.17) and (E.18). After that only one term 
remains: 
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It is the singular part of energy- momentum carried by the electromagnetic field. (It is worth not- 
ing that inverse square root [2E(i, t±, t 2 )\ t2 _ l ._ 00 ]~ 1 / 2 vanishes even if t\ — > — oo, see eqs.(B.ll).) 
Final expressions are presented in Section 4 (see eqs. (4.3.6) and (4.3.7)). 



Appendix G Angular momentum in 2+1 electrodynamics 

We now turn to the calculation of the angular momentum tensor 
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carried by the electromagnetic field due to a pointlike charge. We apply the convenient coordi- 
nate system introduced in Section 4 and detailed in Appendix B. 

We present the integrand of eq.(G.l) in the following form: 
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It is straightforward to substitute the fields (4.1.7) into this expression to calculate the first term 
of the integrand (G.2). The others can be obtained by interchanging of the pair of indices (1, 2) 
and (n, v). 

Having integrated expression Jm^ over ip we obtain 
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where functions A and Ao are given by eqs. (4.2.7). 

Usage of the equalities in eq.(4.2.15) derived in Appendix C allows us to rewrite the integrand 
(G.4) as the following sum: 
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Operators IT 1 are combinations of partial derivatives (4.2.12). 

Further we perform the integration over the time variables and (3 according to the rule in 
eq. (4.3.3). It results in functions of the end points only. We deal with four types of integrals 
described in subsections 4.1°-4.4°. All of them possess a specific small parameter. Near the 
observation time t the small parameter is (5 as well as when t 2 — > — oo. If t 2 tends to h (or 
vice versa), their difference h — t 2 tends to zero. These circumstances simplify the computation 
of integrals of types 4.1°-4.3° which is virtually identical to that presented in Appendix D 
and Appendix E, and we shall not bother with the details. We obtain the bound terms only 
which should be absorbed within the renormalization procedure. Radiative terms arise from the 
integration near the point (3 = (3q where radial variable R = (see integral type 4.4°). 

So, having computed the radiative angular momentum we are not going beyond the limit 
R — ► 0. The terms involved in the final expression (G.5) get simplified sufficiently: 
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They are supplemented with expressions (C.23), (C.38), and (C.34) taken in the point where 
radial variable R = 0. 

First we put \i = and v = % into eq.(G.5). The other terms which constitute the mixed 
spacetime components A^em are obtained by interchanging of indices. A direct consequence of 
the reciprocity is the following combination of partial derivatives in t\ and t 2 : 
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Now we turn to the integration over times t\ and t 2 . It is sufficient to examine the integrals 
near the point R = 0. The computation is virtually identical to that presented in Appendix F. 
After a tedious calculation we obtain the following cumbersome expression: 
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The single integrals belong to the boundary part of angular momentum carried by the elec- 
tromagnetic field. We couple them with integrals over (3 taken at the end points t\ = t 2 and 
t 2 — > — oo. (Such integrals are described in subsections 4.1°-4.4°.) The result is as follows: 
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The double integral in eq.(G.lO) describes the radiative part; it can be rewritten as follows: 
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Taking t 2 — ^ ti limit reveals the proper short-distance behaviour. 
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Now we calculate the space component Mem- Setting fj, = i and v = j into eq.(G.5) we obtain 
Mi2- Having interchanged upper and lower indices we find all the terms which constitute the 
expression .Mem obtained from eq.(G.2) via integration over (p. Further we integrate them over 
times t\ and t 2 . After tedious calculations we finally obtain: 
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All the single integrals should be added to the integrals over (3 evaluated at limit points h = t 2 
and t 2 — ► — oo; the sum is the singular part of angular momentum of the electromagnetic field: 
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The integrand of the radiative part is symmetric in indices (12) and antisymmetric in indices 
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The resulting expressions (G.12) and (G.15) can be rewritten in a manifestly covariant fashion. 
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